3.468 \(\int \frac {(a+b x^3)^{3/2} (c+d x+e x^2+f x^3+g x^4)}{x^7} \, dx\)

Optimal. Leaf size=692 \[ \frac {9\ 3^{3/4} \sqrt {2+\sqrt {3}} b^{2/3} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \left (-5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} b^{2/3} e+4 a g+2 b d\right ) F\left (\sin ^{-1}\left (\frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt [3]{b} x+\left (1+\sqrt {3}\right ) \sqrt [3]{a}}\right )|-7-4 \sqrt {3}\right )}{40 \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}-\frac {27 \sqrt [4]{3} \sqrt {2-\sqrt {3}} \sqrt [3]{a} b^{4/3} e \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} E\left (\sin ^{-1}\left (\frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt [3]{b} x+\left (1+\sqrt {3}\right ) \sqrt [3]{a}}\right )|-7-4 \sqrt {3}\right )}{16 \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}+\frac {27 b^{4/3} e \sqrt {a+b x^3}}{8 \left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )}-\frac {b \sqrt {a+b x^3} \left (10 c x+36 d x^2-45 e x^3-20 f x^4-18 g x^5\right )}{20 x^4}-\frac {1}{60} \left (a+b x^3\right )^{3/2} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )-\frac {b (4 a f+b c) \tanh ^{-1}\left (\frac {\sqrt {a+b x^3}}{\sqrt {a}}\right )}{4 \sqrt {a}}+\frac {b c \sqrt {a+b x^3}}{4 x^3}+\frac {27 b d \sqrt {a+b x^3}}{20 x^2}-\frac {27 b e \sqrt {a+b x^3}}{8 x} \]

[Out]

-1/60*(10*c/x^6+12*d/x^5+15*e/x^4+20*f/x^3+30*g/x^2)*(b*x^3+a)^(3/2)-1/4*b*(4*a*f+b*c)*arctanh((b*x^3+a)^(1/2)
/a^(1/2))/a^(1/2)+1/4*b*c*(b*x^3+a)^(1/2)/x^3+27/20*b*d*(b*x^3+a)^(1/2)/x^2-27/8*b*e*(b*x^3+a)^(1/2)/x-1/20*b*
(-18*g*x^5-20*f*x^4-45*e*x^3+36*d*x^2+10*c*x)*(b*x^3+a)^(1/2)/x^4+27/8*b^(4/3)*e*(b*x^3+a)^(1/2)/(b^(1/3)*x+a^
(1/3)*(1+3^(1/2)))-27/16*3^(1/4)*a^(1/3)*b^(4/3)*e*(a^(1/3)+b^(1/3)*x)*EllipticE((b^(1/3)*x+a^(1/3)*(1-3^(1/2)
))/(b^(1/3)*x+a^(1/3)*(1+3^(1/2))),I*3^(1/2)+2*I)*(1/2*6^(1/2)-1/2*2^(1/2))*((a^(2/3)-a^(1/3)*b^(1/3)*x+b^(2/3
)*x^2)/(b^(1/3)*x+a^(1/3)*(1+3^(1/2)))^2)^(1/2)/(b*x^3+a)^(1/2)/(a^(1/3)*(a^(1/3)+b^(1/3)*x)/(b^(1/3)*x+a^(1/3
)*(1+3^(1/2)))^2)^(1/2)+9/40*3^(3/4)*b^(2/3)*(a^(1/3)+b^(1/3)*x)*EllipticF((b^(1/3)*x+a^(1/3)*(1-3^(1/2)))/(b^
(1/3)*x+a^(1/3)*(1+3^(1/2))),I*3^(1/2)+2*I)*(2*b*d+4*a*g-5*a^(1/3)*b^(2/3)*e*(1-3^(1/2)))*(1/2*6^(1/2)+1/2*2^(
1/2))*((a^(2/3)-a^(1/3)*b^(1/3)*x+b^(2/3)*x^2)/(b^(1/3)*x+a^(1/3)*(1+3^(1/2)))^2)^(1/2)/(b*x^3+a)^(1/2)/(a^(1/
3)*(a^(1/3)+b^(1/3)*x)/(b^(1/3)*x+a^(1/3)*(1+3^(1/2)))^2)^(1/2)

________________________________________________________________________________________

Rubi [A]  time = 1.00, antiderivative size = 692, normalized size of antiderivative = 1.00, number of steps used = 12, number of rules used = 11, integrand size = 35, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.314, Rules used = {14, 1825, 1826, 1835, 1832, 266, 63, 208, 1878, 218, 1877} \[ \frac {9\ 3^{3/4} \sqrt {2+\sqrt {3}} b^{2/3} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \left (-5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} b^{2/3} e+4 a g+2 b d\right ) F\left (\sin ^{-1}\left (\frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt [3]{b} x+\left (1+\sqrt {3}\right ) \sqrt [3]{a}}\right )|-7-4 \sqrt {3}\right )}{40 \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}-\frac {27 \sqrt [4]{3} \sqrt {2-\sqrt {3}} \sqrt [3]{a} b^{4/3} e \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} E\left (\sin ^{-1}\left (\frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt [3]{b} x+\left (1+\sqrt {3}\right ) \sqrt [3]{a}}\right )|-7-4 \sqrt {3}\right )}{16 \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}+\frac {27 b^{4/3} e \sqrt {a+b x^3}}{8 \left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )}-\frac {1}{60} \left (a+b x^3\right )^{3/2} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )-\frac {b \sqrt {a+b x^3} \left (10 c x+36 d x^2-45 e x^3-20 f x^4-18 g x^5\right )}{20 x^4}-\frac {b (4 a f+b c) \tanh ^{-1}\left (\frac {\sqrt {a+b x^3}}{\sqrt {a}}\right )}{4 \sqrt {a}}+\frac {b c \sqrt {a+b x^3}}{4 x^3}+\frac {27 b d \sqrt {a+b x^3}}{20 x^2}-\frac {27 b e \sqrt {a+b x^3}}{8 x} \]

Antiderivative was successfully verified.

[In]

Int[((a + b*x^3)^(3/2)*(c + d*x + e*x^2 + f*x^3 + g*x^4))/x^7,x]

[Out]

(b*c*Sqrt[a + b*x^3])/(4*x^3) + (27*b*d*Sqrt[a + b*x^3])/(20*x^2) - (27*b*e*Sqrt[a + b*x^3])/(8*x) + (27*b^(4/
3)*e*Sqrt[a + b*x^3])/(8*((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)) - (((10*c)/x^6 + (12*d)/x^5 + (15*e)/x^4 + (20*f
)/x^3 + (30*g)/x^2)*(a + b*x^3)^(3/2))/60 - (b*Sqrt[a + b*x^3]*(10*c*x + 36*d*x^2 - 45*e*x^3 - 20*f*x^4 - 18*g
*x^5))/(20*x^4) - (b*(b*c + 4*a*f)*ArcTanh[Sqrt[a + b*x^3]/Sqrt[a]])/(4*Sqrt[a]) - (27*3^(1/4)*Sqrt[2 - Sqrt[3
]]*a^(1/3)*b^(4/3)*e*(a^(1/3) + b^(1/3)*x)*Sqrt[(a^(2/3) - a^(1/3)*b^(1/3)*x + b^(2/3)*x^2)/((1 + Sqrt[3])*a^(
1/3) + b^(1/3)*x)^2]*EllipticE[ArcSin[((1 - Sqrt[3])*a^(1/3) + b^(1/3)*x)/((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)]
, -7 - 4*Sqrt[3]])/(16*Sqrt[(a^(1/3)*(a^(1/3) + b^(1/3)*x))/((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)^2]*Sqrt[a + b*
x^3]) + (9*3^(3/4)*Sqrt[2 + Sqrt[3]]*b^(2/3)*(2*b*d - 5*(1 - Sqrt[3])*a^(1/3)*b^(2/3)*e + 4*a*g)*(a^(1/3) + b^
(1/3)*x)*Sqrt[(a^(2/3) - a^(1/3)*b^(1/3)*x + b^(2/3)*x^2)/((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)^2]*EllipticF[Arc
Sin[((1 - Sqrt[3])*a^(1/3) + b^(1/3)*x)/((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)], -7 - 4*Sqrt[3]])/(40*Sqrt[(a^(1/
3)*(a^(1/3) + b^(1/3)*x))/((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)^2]*Sqrt[a + b*x^3])

Rule 14

Int[(u_)*((c_.)*(x_))^(m_.), x_Symbol] :> Int[ExpandIntegrand[(c*x)^m*u, x], x] /; FreeQ[{c, m}, x] && SumQ[u]
&&  !LinearQ[u, x] &&  !MatchQ[u, (a_) + (b_.)*(v_) /; FreeQ[{a, b}, x] && InverseFunctionQ[v]]

Rule 63

Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> With[{p = Denominator[m]}, Dist[p/b, Sub
st[Int[x^(p*(m + 1) - 1)*(c - (a*d)/b + (d*x^p)/b)^n, x], x, (a + b*x)^(1/p)], x]] /; FreeQ[{a, b, c, d}, x] &
& NeQ[b*c - a*d, 0] && LtQ[-1, m, 0] && LeQ[-1, n, 0] && LeQ[Denominator[n], Denominator[m]] && IntLinearQ[a,
b, c, d, m, n, x]

Rule 208

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-(a/b), 2]*ArcTanh[x/Rt[-(a/b), 2]])/a, x] /; FreeQ[{a,
b}, x] && NegQ[a/b]

Rule 218

Int[1/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = Numer[Rt[b/a, 3]], s = Denom[Rt[b/a, 3]]}, Simp[(2*Sqr
t[2 + Sqrt[3]]*(s + r*x)*Sqrt[(s^2 - r*s*x + r^2*x^2)/((1 + Sqrt[3])*s + r*x)^2]*EllipticF[ArcSin[((1 - Sqrt[3
])*s + r*x)/((1 + Sqrt[3])*s + r*x)], -7 - 4*Sqrt[3]])/(3^(1/4)*r*Sqrt[a + b*x^3]*Sqrt[(s*(s + r*x))/((1 + Sqr
t[3])*s + r*x)^2]), x]] /; FreeQ[{a, b}, x] && PosQ[a]

Rule 266

Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Dist[1/n, Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(a
+ b*x)^p, x], x, x^n], x] /; FreeQ[{a, b, m, n, p}, x] && IntegerQ[Simplify[(m + 1)/n]]

Rule 1825

Int[(Pq_)*(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_.))^(p_), x_Symbol] :> Module[{u = IntHide[x^m*Pq, x]}, Simp[u*(a +
b*x^n)^p, x] - Dist[b*n*p, Int[x^(m + n)*(a + b*x^n)^(p - 1)*ExpandToSum[u/x^(m + 1), x], x], x]] /; FreeQ[{a
, b}, x] && PolyQ[Pq, x] && IGtQ[n, 0] && GtQ[p, 0] && LtQ[m + Expon[Pq, x] + 1, 0]

Rule 1826

Int[(Pq_)*((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_.))^(p_), x_Symbol] :> Module[{q = Expon[Pq, x], i}, Simp[(
c*x)^m*(a + b*x^n)^p*Sum[(Coeff[Pq, x, i]*x^(i + 1))/(m + n*p + i + 1), {i, 0, q}], x] + Dist[a*n*p, Int[(c*x)
^m*(a + b*x^n)^(p - 1)*Sum[(Coeff[Pq, x, i]*x^i)/(m + n*p + i + 1), {i, 0, q}], x], x]] /; FreeQ[{a, b, c, m},
x] && PolyQ[Pq, x] && IGtQ[(n - 1)/2, 0] && GtQ[p, 0]

Rule 1832

Int[(Pq_)/((x_)*Sqrt[(a_) + (b_.)*(x_)^(n_)]), x_Symbol] :> Dist[Coeff[Pq, x, 0], Int[1/(x*Sqrt[a + b*x^n]), x
], x] + Int[ExpandToSum[(Pq - Coeff[Pq, x, 0])/x, x]/Sqrt[a + b*x^n], x] /; FreeQ[{a, b}, x] && PolyQ[Pq, x] &
& IGtQ[n, 0] && NeQ[Coeff[Pq, x, 0], 0]

Rule 1835

Int[(Pq_)*((c_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> With[{Pq0 = Coeff[Pq, x, 0]}, Simp[(Pq
0*(c*x)^(m + 1)*(a + b*x^n)^(p + 1))/(a*c*(m + 1)), x] + Dist[1/(2*a*c*(m + 1)), Int[(c*x)^(m + 1)*ExpandToSum
[(2*a*(m + 1)*(Pq - Pq0))/x - 2*b*Pq0*(m + n*(p + 1) + 1)*x^(n - 1), x]*(a + b*x^n)^p, x], x] /; NeQ[Pq0, 0]]
/; FreeQ[{a, b, c, p}, x] && PolyQ[Pq, x] && IGtQ[n, 0] && LtQ[m, -1] && LeQ[n - 1, Expon[Pq, x]]

Rule 1877

Int[((c_) + (d_.)*(x_))/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = Numer[Simplify[((1 - Sqrt[3])*d)/c]]
, s = Denom[Simplify[((1 - Sqrt[3])*d)/c]]}, Simp[(2*d*s^3*Sqrt[a + b*x^3])/(a*r^2*((1 + Sqrt[3])*s + r*x)), x
] - Simp[(3^(1/4)*Sqrt[2 - Sqrt[3]]*d*s*(s + r*x)*Sqrt[(s^2 - r*s*x + r^2*x^2)/((1 + Sqrt[3])*s + r*x)^2]*Elli
pticE[ArcSin[((1 - Sqrt[3])*s + r*x)/((1 + Sqrt[3])*s + r*x)], -7 - 4*Sqrt[3]])/(r^2*Sqrt[a + b*x^3]*Sqrt[(s*(
s + r*x))/((1 + Sqrt[3])*s + r*x)^2]), x]] /; FreeQ[{a, b, c, d}, x] && PosQ[a] && EqQ[b*c^3 - 2*(5 - 3*Sqrt[3
])*a*d^3, 0]

Rule 1878

Int[((c_) + (d_.)*(x_))/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = Numer[Rt[b/a, 3]], s = Denom[Rt[b/a,
3]]}, Dist[(c*r - (1 - Sqrt[3])*d*s)/r, Int[1/Sqrt[a + b*x^3], x], x] + Dist[d/r, Int[((1 - Sqrt[3])*s + r*x)
/Sqrt[a + b*x^3], x], x]] /; FreeQ[{a, b, c, d}, x] && PosQ[a] && NeQ[b*c^3 - 2*(5 - 3*Sqrt[3])*a*d^3, 0]

Rubi steps

\begin {align*} \int \frac {\left (a+b x^3\right )^{3/2} \left (c+d x+e x^2+f x^3+g x^4\right )}{x^7} \, dx &=-\frac {1}{60} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right ) \left (a+b x^3\right )^{3/2}-\frac {1}{2} (9 b) \int \frac {\sqrt {a+b x^3} \left (-\frac {c}{6}-\frac {d x}{5}-\frac {e x^2}{4}-\frac {f x^3}{3}-\frac {g x^4}{2}\right )}{x^4} \, dx\\ &=-\frac {1}{60} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right ) \left (a+b x^3\right )^{3/2}-\frac {b \sqrt {a+b x^3} \left (10 c x+36 d x^2-45 e x^3-20 f x^4-18 g x^5\right )}{20 x^4}-\frac {1}{4} (27 a b) \int \frac {\frac {c}{9}+\frac {2 d x}{5}-\frac {e x^2}{2}-\frac {2 f x^3}{9}-\frac {g x^4}{5}}{x^4 \sqrt {a+b x^3}} \, dx\\ &=\frac {b c \sqrt {a+b x^3}}{4 x^3}-\frac {1}{60} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right ) \left (a+b x^3\right )^{3/2}-\frac {b \sqrt {a+b x^3} \left (10 c x+36 d x^2-45 e x^3-20 f x^4-18 g x^5\right )}{20 x^4}+\frac {1}{8} (9 b) \int \frac {-\frac {12 a d}{5}+3 a e x+\frac {1}{3} (b c+4 a f) x^2+\frac {6}{5} a g x^3}{x^3 \sqrt {a+b x^3}} \, dx\\ &=\frac {b c \sqrt {a+b x^3}}{4 x^3}+\frac {27 b d \sqrt {a+b x^3}}{20 x^2}-\frac {1}{60} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right ) \left (a+b x^3\right )^{3/2}-\frac {b \sqrt {a+b x^3} \left (10 c x+36 d x^2-45 e x^3-20 f x^4-18 g x^5\right )}{20 x^4}-\frac {(9 b) \int \frac {-12 a^2 e-\frac {4}{3} a (b c+4 a f) x-\frac {12}{5} a (b d+2 a g) x^2}{x^2 \sqrt {a+b x^3}} \, dx}{32 a}\\ &=\frac {b c \sqrt {a+b x^3}}{4 x^3}+\frac {27 b d \sqrt {a+b x^3}}{20 x^2}-\frac {27 b e \sqrt {a+b x^3}}{8 x}-\frac {1}{60} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right ) \left (a+b x^3\right )^{3/2}-\frac {b \sqrt {a+b x^3} \left (10 c x+36 d x^2-45 e x^3-20 f x^4-18 g x^5\right )}{20 x^4}+\frac {(9 b) \int \frac {\frac {8}{3} a^2 (b c+4 a f)+\frac {24}{5} a^2 (b d+2 a g) x+12 a^2 b e x^2}{x \sqrt {a+b x^3}} \, dx}{64 a^2}\\ &=\frac {b c \sqrt {a+b x^3}}{4 x^3}+\frac {27 b d \sqrt {a+b x^3}}{20 x^2}-\frac {27 b e \sqrt {a+b x^3}}{8 x}-\frac {1}{60} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right ) \left (a+b x^3\right )^{3/2}-\frac {b \sqrt {a+b x^3} \left (10 c x+36 d x^2-45 e x^3-20 f x^4-18 g x^5\right )}{20 x^4}+\frac {(9 b) \int \frac {\frac {24}{5} a^2 (b d+2 a g)+12 a^2 b e x}{\sqrt {a+b x^3}} \, dx}{64 a^2}+\frac {1}{8} (3 b (b c+4 a f)) \int \frac {1}{x \sqrt {a+b x^3}} \, dx\\ &=\frac {b c \sqrt {a+b x^3}}{4 x^3}+\frac {27 b d \sqrt {a+b x^3}}{20 x^2}-\frac {27 b e \sqrt {a+b x^3}}{8 x}-\frac {1}{60} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right ) \left (a+b x^3\right )^{3/2}-\frac {b \sqrt {a+b x^3} \left (10 c x+36 d x^2-45 e x^3-20 f x^4-18 g x^5\right )}{20 x^4}+\frac {1}{16} \left (27 b^{5/3} e\right ) \int \frac {\left (1-\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}{\sqrt {a+b x^3}} \, dx+\frac {1}{8} (b (b c+4 a f)) \operatorname {Subst}\left (\int \frac {1}{x \sqrt {a+b x}} \, dx,x,x^3\right )+\frac {1}{80} \left (27 b \left (2 b d-5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} b^{2/3} e+4 a g\right )\right ) \int \frac {1}{\sqrt {a+b x^3}} \, dx\\ &=\frac {b c \sqrt {a+b x^3}}{4 x^3}+\frac {27 b d \sqrt {a+b x^3}}{20 x^2}-\frac {27 b e \sqrt {a+b x^3}}{8 x}+\frac {27 b^{4/3} e \sqrt {a+b x^3}}{8 \left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )}-\frac {1}{60} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right ) \left (a+b x^3\right )^{3/2}-\frac {b \sqrt {a+b x^3} \left (10 c x+36 d x^2-45 e x^3-20 f x^4-18 g x^5\right )}{20 x^4}-\frac {27 \sqrt [4]{3} \sqrt {2-\sqrt {3}} \sqrt [3]{a} b^{4/3} e \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} E\left (\sin ^{-1}\left (\frac {\left (1-\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}{\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}\right )|-7-4 \sqrt {3}\right )}{16 \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}+\frac {9\ 3^{3/4} \sqrt {2+\sqrt {3}} b^{2/3} \left (2 b d-5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} b^{2/3} e+4 a g\right ) \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} F\left (\sin ^{-1}\left (\frac {\left (1-\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}{\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}\right )|-7-4 \sqrt {3}\right )}{40 \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}+\frac {1}{4} (b c+4 a f) \operatorname {Subst}\left (\int \frac {1}{-\frac {a}{b}+\frac {x^2}{b}} \, dx,x,\sqrt {a+b x^3}\right )\\ &=\frac {b c \sqrt {a+b x^3}}{4 x^3}+\frac {27 b d \sqrt {a+b x^3}}{20 x^2}-\frac {27 b e \sqrt {a+b x^3}}{8 x}+\frac {27 b^{4/3} e \sqrt {a+b x^3}}{8 \left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )}-\frac {1}{60} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right ) \left (a+b x^3\right )^{3/2}-\frac {b \sqrt {a+b x^3} \left (10 c x+36 d x^2-45 e x^3-20 f x^4-18 g x^5\right )}{20 x^4}-\frac {b (b c+4 a f) \tanh ^{-1}\left (\frac {\sqrt {a+b x^3}}{\sqrt {a}}\right )}{4 \sqrt {a}}-\frac {27 \sqrt [4]{3} \sqrt {2-\sqrt {3}} \sqrt [3]{a} b^{4/3} e \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} E\left (\sin ^{-1}\left (\frac {\left (1-\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}{\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}\right )|-7-4 \sqrt {3}\right )}{16 \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}+\frac {9\ 3^{3/4} \sqrt {2+\sqrt {3}} b^{2/3} \left (2 b d-5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} b^{2/3} e+4 a g\right ) \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} F\left (\sin ^{-1}\left (\frac {\left (1-\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}{\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}\right )|-7-4 \sqrt {3}\right )}{40 \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}\\ \end {align*}

________________________________________________________________________________________

Mathematica [C]  time = 0.73, size = 240, normalized size = 0.35 \[ \frac {-\frac {12 a^2 d \sqrt {\frac {b x^3}{a}+1} \, _2F_1\left (-\frac {5}{3},-\frac {3}{2};-\frac {2}{3};-\frac {b x^3}{a}\right )}{x^5}-\frac {15 a^2 e \sqrt {\frac {b x^3}{a}+1} \, _2F_1\left (-\frac {3}{2},-\frac {4}{3};-\frac {1}{3};-\frac {b x^3}{a}\right )}{x^4}+\frac {8 b f \left (a+b x^3\right )^3 \, _2F_1\left (2,\frac {5}{2};\frac {7}{2};\frac {b x^3}{a}+1\right )}{a^2}-\frac {30 a^2 g \sqrt {\frac {b x^3}{a}+1} \, _2F_1\left (-\frac {3}{2},-\frac {2}{3};\frac {1}{3};-\frac {b x^3}{a}\right )}{x^2}-15 b^2 c \sqrt {\frac {b x^3}{a}+1} \tanh ^{-1}\left (\sqrt {\frac {b x^3}{a}+1}\right )-\frac {15 b c \left (a+b x^3\right )}{x^3}-\frac {10 c \left (a+b x^3\right )^2}{x^6}}{60 \sqrt {a+b x^3}} \]

Antiderivative was successfully verified.

[In]

Integrate[((a + b*x^3)^(3/2)*(c + d*x + e*x^2 + f*x^3 + g*x^4))/x^7,x]

[Out]

((-15*b*c*(a + b*x^3))/x^3 - (10*c*(a + b*x^3)^2)/x^6 - 15*b^2*c*Sqrt[1 + (b*x^3)/a]*ArcTanh[Sqrt[1 + (b*x^3)/
a]] - (12*a^2*d*Sqrt[1 + (b*x^3)/a]*Hypergeometric2F1[-5/3, -3/2, -2/3, -((b*x^3)/a)])/x^5 - (15*a^2*e*Sqrt[1
+ (b*x^3)/a]*Hypergeometric2F1[-3/2, -4/3, -1/3, -((b*x^3)/a)])/x^4 - (30*a^2*g*Sqrt[1 + (b*x^3)/a]*Hypergeome
tric2F1[-3/2, -2/3, 1/3, -((b*x^3)/a)])/x^2 + (8*b*f*(a + b*x^3)^3*Hypergeometric2F1[2, 5/2, 7/2, 1 + (b*x^3)/
a])/a^2)/(60*Sqrt[a + b*x^3])

________________________________________________________________________________________

fricas [F]  time = 0.99, size = 0, normalized size = 0.00 \[ {\rm integral}\left (\frac {{\left (b g x^{7} + b f x^{6} + b e x^{5} + {\left (b d + a g\right )} x^{4} + a e x^{2} + {\left (b c + a f\right )} x^{3} + a d x + a c\right )} \sqrt {b x^{3} + a}}{x^{7}}, x\right ) \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b*x^3+a)^(3/2)*(g*x^4+f*x^3+e*x^2+d*x+c)/x^7,x, algorithm="fricas")

[Out]

integral((b*g*x^7 + b*f*x^6 + b*e*x^5 + (b*d + a*g)*x^4 + a*e*x^2 + (b*c + a*f)*x^3 + a*d*x + a*c)*sqrt(b*x^3
+ a)/x^7, x)

________________________________________________________________________________________

giac [F]  time = 0.00, size = 0, normalized size = 0.00 \[ \int \frac {{\left (g x^{4} + f x^{3} + e x^{2} + d x + c\right )} {\left (b x^{3} + a\right )}^{\frac {3}{2}}}{x^{7}}\,{d x} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b*x^3+a)^(3/2)*(g*x^4+f*x^3+e*x^2+d*x+c)/x^7,x, algorithm="giac")

[Out]

integrate((g*x^4 + f*x^3 + e*x^2 + d*x + c)*(b*x^3 + a)^(3/2)/x^7, x)

________________________________________________________________________________________

maple [B]  time = 0.07, size = 1196, normalized size = 1.73 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((b*x^3+a)^(3/2)*(g*x^4+f*x^3+e*x^2+d*x+c)/x^7,x)

[Out]

d*(-1/5*(b*x^3+a)^(1/2)*a/x^5-13/20*(b*x^3+a)^(1/2)*b/x^2-9/20*I*b*3^(1/2)*(-a*b^2)^(1/3)*(I*(x+1/2*(-a*b^2)^(
1/3)/b-1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b)*3^(1/2)/(-a*b^2)^(1/3)*b)^(1/2)*((x-(-a*b^2)^(1/3)/b)/(-3/2*(-a*b^2)^(1
/3)/b+1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b))^(1/2)*(-I*(x+1/2*(-a*b^2)^(1/3)/b+1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b)*3^(1/
2)/(-a*b^2)^(1/3)*b)^(1/2)/(b*x^3+a)^(1/2)*EllipticF(1/3*3^(1/2)*(I*(x+1/2*(-a*b^2)^(1/3)/b-1/2*I*3^(1/2)*(-a*
b^2)^(1/3)/b)*3^(1/2)/(-a*b^2)^(1/3)*b)^(1/2),(I*3^(1/2)*(-a*b^2)^(1/3)/(-3/2*(-a*b^2)^(1/3)/b+1/2*I*3^(1/2)*(
-a*b^2)^(1/3)/b)/b)^(1/2)))+e*(-1/4*(b*x^3+a)^(1/2)*a/x^4-11/8*(b*x^3+a)^(1/2)*b/x-9/8*I*b*3^(1/2)*(-a*b^2)^(1
/3)*(I*(x+1/2*(-a*b^2)^(1/3)/b-1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b)*3^(1/2)/(-a*b^2)^(1/3)*b)^(1/2)*((x-(-a*b^2)^(1
/3)/b)/(-3/2*(-a*b^2)^(1/3)/b+1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b))^(1/2)*(-I*(x+1/2*(-a*b^2)^(1/3)/b+1/2*I*3^(1/2)
*(-a*b^2)^(1/3)/b)*3^(1/2)/(-a*b^2)^(1/3)*b)^(1/2)/(b*x^3+a)^(1/2)*((-3/2*(-a*b^2)^(1/3)/b+1/2*I*3^(1/2)*(-a*b
^2)^(1/3)/b)*EllipticE(1/3*3^(1/2)*(I*(x+1/2*(-a*b^2)^(1/3)/b-1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b)*3^(1/2)/(-a*b^2)
^(1/3)*b)^(1/2),(I*3^(1/2)*(-a*b^2)^(1/3)/(-3/2*(-a*b^2)^(1/3)/b+1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b)/b)^(1/2))+(-a
*b^2)^(1/3)/b*EllipticF(1/3*3^(1/2)*(I*(x+1/2*(-a*b^2)^(1/3)/b-1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b)*3^(1/2)/(-a*b^2
)^(1/3)*b)^(1/2),(I*3^(1/2)*(-a*b^2)^(1/3)/(-3/2*(-a*b^2)^(1/3)/b+1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b)/b)^(1/2))))+
f*(-1/3*(b*x^3+a)^(1/2)*a/x^3+2/3*(b*x^3+a)^(1/2)*b-b*arctanh((b*x^3+a)^(1/2)/a^(1/2))*a^(1/2))+g*(-1/2*(b*x^3
+a)^(1/2)*a/x^2+2/5*(b*x^3+a)^(1/2)*b*x-9/10*I*a*3^(1/2)*(-a*b^2)^(1/3)*(I*(x+1/2*(-a*b^2)^(1/3)/b-1/2*I*3^(1/
2)*(-a*b^2)^(1/3)/b)*3^(1/2)/(-a*b^2)^(1/3)*b)^(1/2)*((x-(-a*b^2)^(1/3)/b)/(-3/2*(-a*b^2)^(1/3)/b+1/2*I*3^(1/2
)*(-a*b^2)^(1/3)/b))^(1/2)*(-I*(x+1/2*(-a*b^2)^(1/3)/b+1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b)*3^(1/2)/(-a*b^2)^(1/3)*
b)^(1/2)/(b*x^3+a)^(1/2)*EllipticF(1/3*3^(1/2)*(I*(x+1/2*(-a*b^2)^(1/3)/b-1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b)*3^(1
/2)/(-a*b^2)^(1/3)*b)^(1/2),(I*3^(1/2)*(-a*b^2)^(1/3)/(-3/2*(-a*b^2)^(1/3)/b+1/2*I*3^(1/2)*(-a*b^2)^(1/3)/b)/b
)^(1/2)))+c*(-1/6*a*(b*x^3+a)^(1/2)/x^6-5/12*(b*x^3+a)^(1/2)*b/x^3-1/4*b^2*arctanh((b*x^3+a)^(1/2)/a^(1/2))/a^
(1/2))

________________________________________________________________________________________

maxima [F]  time = 0.00, size = 0, normalized size = 0.00 \[ \frac {1}{24} \, {\left (\frac {3 \, b^{2} \log \left (\frac {\sqrt {b x^{3} + a} - \sqrt {a}}{\sqrt {b x^{3} + a} + \sqrt {a}}\right )}{\sqrt {a}} - \frac {2 \, {\left (5 \, {\left (b x^{3} + a\right )}^{\frac {3}{2}} b^{2} - 3 \, \sqrt {b x^{3} + a} a b^{2}\right )}}{{\left (b x^{3} + a\right )}^{2} - 2 \, {\left (b x^{3} + a\right )} a + a^{2}}\right )} c + \int \frac {{\left (b g x^{6} + b f x^{5} + b e x^{4} + a f x^{2} + {\left (b d + a g\right )} x^{3} + a e x + a d\right )} \sqrt {b x^{3} + a}}{x^{6}}\,{d x} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b*x^3+a)^(3/2)*(g*x^4+f*x^3+e*x^2+d*x+c)/x^7,x, algorithm="maxima")

[Out]

1/24*(3*b^2*log((sqrt(b*x^3 + a) - sqrt(a))/(sqrt(b*x^3 + a) + sqrt(a)))/sqrt(a) - 2*(5*(b*x^3 + a)^(3/2)*b^2
- 3*sqrt(b*x^3 + a)*a*b^2)/((b*x^3 + a)^2 - 2*(b*x^3 + a)*a + a^2))*c + integrate((b*g*x^6 + b*f*x^5 + b*e*x^4
+ a*f*x^2 + (b*d + a*g)*x^3 + a*e*x + a*d)*sqrt(b*x^3 + a)/x^6, x)

________________________________________________________________________________________

mupad [F]  time = 0.00, size = -1, normalized size = -0.00 \[ \int \frac {{\left (b\,x^3+a\right )}^{3/2}\,\left (g\,x^4+f\,x^3+e\,x^2+d\,x+c\right )}{x^7} \,d x \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(((a + b*x^3)^(3/2)*(c + d*x + e*x^2 + f*x^3 + g*x^4))/x^7,x)

[Out]

int(((a + b*x^3)^(3/2)*(c + d*x + e*x^2 + f*x^3 + g*x^4))/x^7, x)

________________________________________________________________________________________

sympy [A]  time = 17.92, size = 524, normalized size = 0.76 \[ \frac {a^{\frac {3}{2}} d \Gamma \left (- \frac {5}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {5}{3}, - \frac {1}{2} \\ - \frac {2}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 x^{5} \Gamma \left (- \frac {2}{3}\right )} + \frac {a^{\frac {3}{2}} e \Gamma \left (- \frac {4}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {4}{3}, - \frac {1}{2} \\ - \frac {1}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 x^{4} \Gamma \left (- \frac {1}{3}\right )} + \frac {a^{\frac {3}{2}} g \Gamma \left (- \frac {2}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {2}{3}, - \frac {1}{2} \\ \frac {1}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 x^{2} \Gamma \left (\frac {1}{3}\right )} + \frac {\sqrt {a} b d \Gamma \left (- \frac {2}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {2}{3}, - \frac {1}{2} \\ \frac {1}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 x^{2} \Gamma \left (\frac {1}{3}\right )} + \frac {\sqrt {a} b e \Gamma \left (- \frac {1}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {1}{2}, - \frac {1}{3} \\ \frac {2}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 x \Gamma \left (\frac {2}{3}\right )} - \sqrt {a} b f \operatorname {asinh}{\left (\frac {\sqrt {a}}{\sqrt {b} x^{\frac {3}{2}}} \right )} + \frac {\sqrt {a} b g x \Gamma \left (\frac {1}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {1}{2}, \frac {1}{3} \\ \frac {4}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 \Gamma \left (\frac {4}{3}\right )} - \frac {a^{2} c}{6 \sqrt {b} x^{\frac {15}{2}} \sqrt {\frac {a}{b x^{3}} + 1}} - \frac {a \sqrt {b} c}{4 x^{\frac {9}{2}} \sqrt {\frac {a}{b x^{3}} + 1}} - \frac {a \sqrt {b} f \sqrt {\frac {a}{b x^{3}} + 1}}{3 x^{\frac {3}{2}}} + \frac {2 a \sqrt {b} f}{3 x^{\frac {3}{2}} \sqrt {\frac {a}{b x^{3}} + 1}} - \frac {b^{\frac {3}{2}} c \sqrt {\frac {a}{b x^{3}} + 1}}{3 x^{\frac {3}{2}}} - \frac {b^{\frac {3}{2}} c}{12 x^{\frac {3}{2}} \sqrt {\frac {a}{b x^{3}} + 1}} + \frac {2 b^{\frac {3}{2}} f x^{\frac {3}{2}}}{3 \sqrt {\frac {a}{b x^{3}} + 1}} - \frac {b^{2} c \operatorname {asinh}{\left (\frac {\sqrt {a}}{\sqrt {b} x^{\frac {3}{2}}} \right )}}{4 \sqrt {a}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((b*x**3+a)**(3/2)*(g*x**4+f*x**3+e*x**2+d*x+c)/x**7,x)

[Out]

a**(3/2)*d*gamma(-5/3)*hyper((-5/3, -1/2), (-2/3,), b*x**3*exp_polar(I*pi)/a)/(3*x**5*gamma(-2/3)) + a**(3/2)*
e*gamma(-4/3)*hyper((-4/3, -1/2), (-1/3,), b*x**3*exp_polar(I*pi)/a)/(3*x**4*gamma(-1/3)) + a**(3/2)*g*gamma(-
2/3)*hyper((-2/3, -1/2), (1/3,), b*x**3*exp_polar(I*pi)/a)/(3*x**2*gamma(1/3)) + sqrt(a)*b*d*gamma(-2/3)*hyper
((-2/3, -1/2), (1/3,), b*x**3*exp_polar(I*pi)/a)/(3*x**2*gamma(1/3)) + sqrt(a)*b*e*gamma(-1/3)*hyper((-1/2, -1
/3), (2/3,), b*x**3*exp_polar(I*pi)/a)/(3*x*gamma(2/3)) - sqrt(a)*b*f*asinh(sqrt(a)/(sqrt(b)*x**(3/2))) + sqrt
(a)*b*g*x*gamma(1/3)*hyper((-1/2, 1/3), (4/3,), b*x**3*exp_polar(I*pi)/a)/(3*gamma(4/3)) - a**2*c/(6*sqrt(b)*x
**(15/2)*sqrt(a/(b*x**3) + 1)) - a*sqrt(b)*c/(4*x**(9/2)*sqrt(a/(b*x**3) + 1)) - a*sqrt(b)*f*sqrt(a/(b*x**3) +
1)/(3*x**(3/2)) + 2*a*sqrt(b)*f/(3*x**(3/2)*sqrt(a/(b*x**3) + 1)) - b**(3/2)*c*sqrt(a/(b*x**3) + 1)/(3*x**(3/
2)) - b**(3/2)*c/(12*x**(3/2)*sqrt(a/(b*x**3) + 1)) + 2*b**(3/2)*f*x**(3/2)/(3*sqrt(a/(b*x**3) + 1)) - b**2*c*
asinh(sqrt(a)/(sqrt(b)*x**(3/2)))/(4*sqrt(a))

________________________________________________________________________________________